DIRECT AND INVERSE OBSTACLE SCATTERING PROBLEMS IN 
A PIECEWISE HOMOGENEOUS MEDIUM 
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Abstract. This paper is concerned with the problem of scattering of time-harmonic acoustic 
waves from an impenetrable obstacle in a piecewise homogeneous medium. The well-posedness of the 
direct problem is established, employing the integral equation method and then used, in conjunction 
with the representation in a combination of layer potentials of the solution, to prove a priori estimates 
of solutions on some part of the interface between the layered media. The inverse problem is also 
considered in this paper. An uniqueness result is obtained for the first time in determining both the 
penetrable interface and the impenetrable obstacle with its physical property from a knowledge of 
the far field pattern for incident plane waves. In doing so, an important role is played by the a priori 
estimates of the solution for the direct problem. 
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1. Introduction. In this paper, we consider the problem of scattering of time- 
harmonic acoustic plane waves by an impenetrable obstacle in a piecewise homoge- 
neous medium. In practical applications, the background might not be homogeneous 
and then may be modeled as a layered medium. A medium of this type that is a 
nested body consisting of a finite number of homogeneous layers occurs in various 
areas of applications such as radar, remote sensing, geophysics, and nondestructive 
testing. 

To give a precise description of the problem, let C K."^ denote the impenetrable 
obstacle which is an open bounded region with a boundary Si and let E'^\02 
denote the the background medium which is divided by means of a closed surface 
So into two connected domains and fli (see Figure [TTT|) . Here, fio is the unbounded 
homogeneous medium and fli is the bounded homogeneous one. We assume that the 
boundary Si of the obstacle fl2 has a dissection Si = Fq U Fi, where Fq and Fi are 
two disjoint, relatively open subsets of Si. Furthermore, the Dirichlet and impedance 
boundary conditions with the surface impedance a nonnegative continuous function 
A € C'(Fi) are specified on Fq and Fi, respectively. Note that the case Fi = 
corresponds to a sound-soft obstacle and the case Fq = 0, A = leads to a Neumann 
boundary condition which corresponds to a sound-hard obstacle. 

The scattering of time-harmonic acoustic waves in a two-layered medium in is 
now modeled by the Helmholtz equation with boundary conditions on the interface 
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Sq and boundary 5*1: 
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where v is the unit outward normal to the interface 5*0 and boundary 5*1, Xq is a 
positive constant. Here, the total field u — +u' is given as the sum of the unknown 
scattered wave u" which is required to satisfy the Sommerfeld radiation condition (|1.5p 
and incident plane wave = e^^ox-d ^ where kj is the positive wave number given by 
kj — ujj/cj in terms of the frequency Wj and the sound speed Cj in the corresponding 
region Vlj (j — 0, 1). The distinct wave numbers kj {j = 0, 1) correspond to the fact 
that the background medium consists of two physically different materials. On the 
interface 5*0, the so-called "transmission condition" (|1.3p is imposed, which represents 
the continuity of the medium and equilibrium of the forces acting on it. The boundary 
condition 3i{v) = on S'o is understood as: 

v = on To, (1.6) 

dv 

— +iXv = onTi. (1.7) 
ov 

Thus, the boundary condition (|1.4p is a general and realistic one and allows that the 
pressure of the total wave v vanishes on Fg and the normal velocity is proportional 
to the excess pressure on the coated part Fi. 




Fig. 1.1. Scattering in a two-layered background medium 

The direct problem is to seek a pair of functions u S C^(r2o) H C"'^'"(Slo) and 
V e C2(0i) n C^''^(Jh) satisfying ((Ll])-(ITl5l). By the variational method, the well- 
posedness (existence, uniqueness and stability) of the direct problem has been es- 
tablished in 'T for the Dirichlet boundary condition and in [18] for a general mixed 
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boundary condition (|1.4p . In the present paper, an integral equation method is em- 
ployed to establish the well-posedness of the direct problem. This result is also used, 
in conjunction with the representation in a combination of layer potentials of the so- 
lution, to prove a priori estimates of the solution on some part of the interface 5*0, 
which plays an important role in the proof of the uniqueness result for our inverse 
problem later on. 

Further, it is known that ^^(x) has the following asymptotic representation 



uniformly for all directions x := x/\x\, where the function u°°(x,d) defined on the 
unit sphere S is known as the far field pattern with x and d denoting, respectively, 
the observation direction and the incident direction. 

The inverse problem we consider in this paper is, given the wave numbers kj 
(j ~ 0, 1), the positive constant Aq and the far field pattern u°°{x, d) for all incident 
plane waves with incident direction d £ S*, to determine the obstacle with its 
physical property and the interface 5*0. As usual in most of the inverse problems, 
the first question to ask in this context is the identifiability, that is, whether an 
inaccessible obstacle 0,2 with its physical property i3§ and the interface 5o can be 
identified from a knowledge of the far-field pattern. Mathematically, the identifiability 
is the uniqueness issue which is of theoretical interest and is required in order to 
proceed to efficient numerical methods of solutions. 

Since the first uniqueness result given by Schiffer in 1967 for a sound-soft obstacle 
[HITS], there has been an extensive study in this direction in the literature; see, e.g. 
[3[3[3[S[T0l[lll[20l[2a|23|2l[25l[26l[^for scattering in a homogeneous medium, 
O [131 [19] for scattering in an inhomogeneous medium and pjj for scattering by special 
obstacles such as balls and polyhedra. However, there are few uniqueness results for 
inverse obstacle scattering in a piecewise homogeneous medium. For the case of a 
known piecewise homogeneous medium, Yan and Pang [29) established a uniqueness 
result for the inverse scattering problem of determining a sound-soft obstacle based on 
Schiffer's idea; their method can not be extended to other boundary conditions. They 
obtained a uniqueness result for the case of a sound-hard obstacle in a two-layered 
background medium in [5T] using a generalization of Schiffer's method. However, their 
method is hard to be extended to the case of a multilayered background medium and 
seems unreasonable to require the interior wave number to be in an interval. Recently 
in [18j . based on a generalization of the mixed reciprocity relation, we proved that 
both the obstacle D,2 and its physical property S§ can be uniquely recovered from 
a knowledge of the far field pattern for incident plane waves. This seems to be 
appropriate for a number of applications where the physical nature of the obstacle is 
unknown. The tools and the uniqueness result developed in [TS] can also be extended 
to inverse electromagnetic scattering problems |17j . For the case of an unknown 
piecewise homogeneous medium, Athanasiadis, Ramm and Stratis 1 and Yan |28| 
proved that the interfaces between the layered media can be determined uniquely by 
the corresponding far field pattern in the special case when the impenetrable obstacle 
does not exist. 

However, to the authors' knowledge, no uniqueness result is available for deter- 
mining both the obstacle embedded in the piecewise homogeneous medium and the 
interfaces between the layered media from a knowledge of the far field pattern for 
incident plane waves. In this paper, we will prove for the first time that both the 
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inaccessible obstacle with its physical property and the interface 5*0 can be 
uniquely determined by a knowledge of the far-field pattern. We remark that the 
results obtained in this paper are also available for both the 2D case and the case of 
a multilayered medium and can be proved similarly. 

The remaining part of the paper is organized as follows. In the next section, 
we will establish the well-posedness of the direct scattering problem, employing the 
integral equation method. With the help of the representation in a combination of 
layer potentials of the solution, a priori estimates of solutions are also obtained on 
some part of the interface between the layered media. Section[3]is devoted to the proof 
of the result on the unique determination of both the obstacle with its physical 
property ^ and the surface Sq from a knowledge of the far field pattern for incident 
plane waves. 

2. The direct scattering problem. In this section we first establish the well- 
posedness of the direct problem, employing the integral equation method and then 
make use of the representation in a combination of layer potentials of the solution 
to prove some a priori estimates of the solution which plays an important role in 
the inverse problem. We shall use C to denote a generic constant whose values may 
change in different inequalities but always bounded away from infinity. 

As incident fields u', plane waves and point sources (cf. (|2.7p below) are of special 
interest. Denote by d) the scattered field for an incident plane wave rf) with 
incident direction d G S and by u°°{-,d) the corresponding far field pattern. The 
scattered field for an incident point source <&(•, with source point z G is denoted 
by u*(-; z) and the corresponding far field pattern by $°°(-, z). 

The direct problem is to look for a pair of functions u e C^(rio) n C^'"(r2o) and 
V e C^(r2i) n C^'"{Qi) satisfying the following boundary value problem: 



Au + klu^O inflo, (2.1) 

Av + klv = Q inf^i, (2.2) 

du dv 

u-v = f, AoT^=.g on So, (2.3) 

oi' ov 

v^Q on To, (2.4) 

-^+i\v = {) onTi, (2.5) 
ov 

du 

lim r(— ikou) = 0,r = |a;|. (2-6) 

r— >oo Of 



Here, we assume that fcp, ki and Aq are given positive constants and that / G C^'"(S'o) 
and g G C°'"(S'o) are given functions in Holder spaces with exponent < a < 1. 

Remark 2.1. The problem of scattering of the incident plane wave = e^^"^'"^ 
is a particular case of the problem (|2.ip - (|2.6p . In particular, the scattered field u" 

satisfies the problem (I2.ip - (l2.6p with , / = — u'k g — —-;:^\sa- 

ov 

The following uniqueness result has been established in [TS] (see Theorem 2.3 of 

m)- 

Theorem 2.2. The boundary value problem ^2.1\ — \2.b\ admits at most one 
solution. 

Denote by $j the fundamental solution of the Helmholtz equation with wave 
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number kj {j — 0,1), which is given by 

^ikj\x-y\ 

'^j{x:y) = ir\ T' x,yeR\xj^y. (2.7) 

47r|a; — y\ 

For i,j = 0, 1 define the single- and double-layer operators Si j and Ki j, respectively, 

by 

{S,^j(j)){x) := ^j{x,y)(f){y)ds{y) x e S,, 



(K,,M^) := / ^^^^cl,{y)ds{y) x € 
Jso c>v{y) 



and the normal derivative operators j and Tij by 



d 



{K,,j(l)){x) -.^ <^j{x,y)(l){y)ds{y) x e S,. 



d f d<^>j{x,y) 



For i,j = 0, 1 define the single- and double-layer operators Si j and Ki j, respectively, 

by 

iS,,j<j>)ix):^ f ^i{x,y)<i>{y)ds(y) x & S,, (2.8) 



iK^,,c^)ix) := ^^^c^{y)ds{y) x 6 S. (2.9) 



and the normal derivative operators K- ■ and Ti^ by 

{kIj4>){x) ^^-j-^ j ^i{x,y)(t>{y)ds{y) x e S,. 



dv{x) 



r. 



{T,,M^) ^ / ^il^0(y)ds(2/) X e 5. 

For the mapping properties of these operators in the spaces of continuous and 
Holder continuous functions we refer to Section 3.1 in [4j or Chapter 2 in [3]. 

Theorem 2.3. The boundary value problem ^2.1\ — \2.6^ has a unique solution. 
Moreover, there exists a positive constant C — C{\Q,U,i,a) such that 

Mc^-^m + lkllci^=(or) ^ c{\\f\\c^^.(So) + ll5llco-(5o))- (2.10) 



Proof. The uniqueness of solutions follows from Theorem 12.21 We now prove the 
existence of solutions by using the integral equation method. Following [3] and [3] we 
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seek a solution in the form 

d<^>o{x,y) 



u{x) 
v{x) 



So 



ro 



d^i{x,y) 
dv{y) 

d^i{x,y) 
du{y) 



i:{y) + ^o{x,y)(t){y)\ ds{y) 



xG^o, (2.11) 



^{y) + <i?iix,y)4:{y) j ds{y) 
X{y) - ir]^i{x,y)x{y) \ ds{y) 



^iix,y)ip{y) +iri 



d<t>i{x,y ) ^ 
dv{y) 



{S''^){y))\ds{y) xeni (2.12) 



with four densities ?A e C^-^iSo), e^C"'"(S'o), x e ^^^"(ro), ip £ C^''^{Ti) and 
a real coupling parameter rj ^ 0. By 5* we denote the single-layer operator (|2.8p in 
the potential theoretic limit case ki = 0. Then from the jump relations we see that 
the potentials u and v defined above solve the boundary value problem (|2.ip — (|2.6p 
provided the densities tp, 0, %, if satisfy the system of integral equations 



-^(7^0,0 - ivSo.o)x + A/i(y5 = nf on 5*0, 
(j) - /iAo(ro.o - T()s)if^ - fi{K„ Q - Aoifoa)0 

~^Xo{Tnfi - i'nKQ Q)x + = -fig on ^o, 
X + 2Xi,iV' + 2S'i,i0 + 2(Xi,o - ivSi,o)x + Msip = on Tq, 



(2.13) 

(2.14) 
(2.15) 

LP - 2(ri,i + iXKi^i)^} - 2(^1 1 + iXSi^i)(j) + AUx + M^Lp = on Ti (2.16) 
with /i = 2/(Ao + 1), where 

Ml = -M(5oa+ir;i^o,i^'), 

M2 = -M(/^oa+i^?o,i^'), 

Ma = 2(^1,1 +Z77i^i,i^'), 

M4 - -2(fi,o - ii^k[ a + aXi,o + A?7^i,o), 

M5 = -2(^14 + ir]fi^iS^ - XrjS'^ + iXSi,i ~ Xr^Ki.iS^). 

Define the product space X := C^^^iSo) x CO^"(S'o) x C^^^iVo) x Ci'"(ri) and 
introduce the operator A : X ^ X given by 



A := 



-mAo(7o,o - To,i) -^(Xq o - AoXo.i) -f^^o(To,o - itiKq q) M2 
2ifia 2^1,1 ' 2(7^1,0 - z77^i^o) M3 

\-2(Ti,i + 2Aifi,i) -2(if;i + 2AS'i,i) M4 M5/ 



The operator ^ is compact since all its entries are compact. The system (|2.13p - (|2.16p 
can be rewritten in the abbreviated form 



iI + A)U^R, 



(2.17) 



where / is the identity operator, U — (tp, (/), Xi V)^ ^^^d = (m/j ~M5i Oj 0)"^- Thus, 
the Riesz-Fredholm theory is applicable. We now prove the uniqueness of solutions 
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to the system (|2.17p . To this end, let [/ be a solution of the homogeneous system 
corresponding to (|2.17p (that is, the system (|2.17p with R ~ 0). Then it is enough to 
show that U = 0. 

We first prove that x = on Fq and = on Fi. From the system (I2.f 7p or 
(|2.13p - (|2.16p with fif ~ —fig = (since i? = 0) it is known that u and v defined 
in ([2TT|l and ([TT2| satisfy the problem - dm with f = g = 0. Thus, by the 
uniqueness Theorem 12.21 u = in fio and u = in fii. Note that v, given by (|2.f 2p . 
can also be defined for x E and satisfies the Helmholtz equation Av + kfv = in 
Q2- Then the jump relations yield that 

dv- . 

-V- = X, — ^ = ^VX on Fq, 
—V- — iriS </J, — = —CO on f 1. 

Of 

Interchanging the order of integration and using Green's first theorem over 0,2, we 
obtain 

/ \x?ds+ / \SLp\^ds\ ^ i-ql / 1x1^^5+ / (pS^^ds 

dv^ f _ dv^ 

— as + / — as 

av Jn ov 



[ \Vv\^ ~ kl\v\^dx. 



Taking the imaginary part of this equation gives that x = on Fq and Sip = on Fi. 
The single-layer potential 

"^^^^ ■= / TT~ -My)ds{y) 

Jr, 47r|x - y\ 

with density (p is continuous throughout MP, harmonic in R'^yFi and vanishes on 
Fi and at infinity. Therefore, by the maximum-minimum principle for harmonic 
functions, we have w = in M'^ and the jump relation yields (p = 0. 
Now the system (|2.I7p becomes 







Define 



v{^)--= l^^^^^j{y) + 'S>i{x,y)4>{y)^ds{y), x e O^, 

u{x):=- |^^^V(2/) + ^<i>o(x,2/)</>(y)|ds(y), x &M?\Tl^. 

Then by the jump relations for single- and double-layer potentials we have 

v — v — ip, —u + u = %jj on 5*0, (2-18) 
Ao 

dv dv ^ du ^ du , „ _ 

7r-7r = -'^' 7r + ^07r = -'^ °" -^n- (2.19) 
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Hence, v and u solve the homogeneous transmission problem 

Av + klv = in fip, Au + klu = in R^\IIo 

with the transmission conditions (noting that u = in fJo and w = in f2i) 

~ ~ n dv du 

v-u = 0, — = Ao-^ on Sq. 

Of Of 

Arguing similarly as in the proof of Theorem 2.3 in |18j we can show that u = in 
^^0 and w = in R'^\ilo. Hence we conclude from (I2.18|) and (|2.19p that ^ = (/) — 
on 5*0. 

Thus, the injectivity of the operator I + A is proved, and by the Riesz-Fredholm 
theory {I + A)~^ exists and is bounded in X. From this we deduce that 

IIV'llci.°(So) + U\\c°-''{So) + llxllci>°(ro) + ll¥'llci.°(ri) 

<C(||/||ci.= (So) + ll.9llco^»(5„))- 

The estimate (|2l0l) follows from ([2JT|) . (12^2]) and Theorem 3.3 in L4J. □ 

We now make use of the representation (I2.12p of the solution v to derive an a 

priori estimate of the solution v on some part of 5*0, which is necessary in proving the 

uniqueness result for the inverse problem in the next section. 

Let X* G 5*0 be an arbitrarily fixed point and let us introduce the space Cq{So) 

which consists of all continuous functions h G C(S'o\{a::*}) with the property that 

lim [\x — x*\h(x)] 

X — ^x* 

exists. It can easily be seen that Co(<S'o) is a Banach space equipped with the weighted 
maximum norm 

II^IICo(So) — sup \{x - X*)h{x)\. 

x^x* ,x^Sq 

Lemma 2.4. Given two functions f S C^'"(5'o) and g e C°'"(S'o). Let u e 
C2(r2o) n C^'°'(Sh) and V G C^iVLi) n Ci'"(Or) be the solution of the problem ^EM - 
(|i?.6'p . Let X* G So and let Bi, B2 be two small balls with center x* and radii ri, r2, 
respectively, satisfying that ri < 7'2 . Then there exists a constant C > such that 



M 



',So\B2 



dv 
dv 



oo,So\-B2 

< ^^(11/1100(50) + ll5llco(So) + ll/lll,a,So\Si + ll.9llo,a,5o\Bi ) 



Proof. We consider again the system (|2.17p of the boundary integral equations 
derived from the boundary value problem (|2.ip - p.6p . In addition to the space X, 
we also consider the weighted spaces Co Co(5o) x Co(S'o) x C^'"(ro) x C^'"(ri). 
The matrix operator A is also compact in Co since all entries of A are compact (see 
[I2I [13]). From the proof of Theorem 12. 3i we know that the operator I + A has a 
trivial null space in X. Therefore, by the Fredholm alternative applied to the dual 
system (X, Co) with the L^ bilinear form, the adjoint operator I + A' has a trivial null 
space in Co. By the Fredholm alternative again, but now applied to the dual system 
(Co, Co) with the L^ bilinear form, the operator I + A also has a trivial null space in 
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Cq. Hence, by the Riesz-Fredholm theory, the system (|2.17l) is also uniquely solvable 
in Co, and the solution depends continuously on the right-hand side: 

IIV'llco(So) + ll<^llco(So) + llxllci>=(ro) + ll'Pllci.°(ri) 

< C(||/||co(So) + ll3llco(So))- (2.20) 

From (|2.12p and the jump relation we find that on Sq, 

V = -i^A + Ko,i^ + Sq.kI) + {kofl - iilSo.o)x + {Sq,i + ivKQsS^)ip. (2.21) 
Thus we have 

lklloo,So\i32 < C'(IIV'Iloo,So\B2 + ll^''o4V'lloo,So\B2 + l|5'o4</'Iloo,So\B2 

+ 11(^0,0 - ■ivSofi)x\\oo,So\B2 + \\{SoS + ■ivKQ,lS^)(p\\oo,So\B2) 
Co{So) ~r" 11^0,1^11 

+ llxllci-(ro) + ll^llci.»(ri))- (2-22) 

We choose a function pi e C^(S'o) such that pi{x) = for x G Sq\Bi and pi{x) = 1 
in the neighborhood of x* . We also choose another function p2 £ C^(5o) such that 
P2{x) = 1 for a; G S'o\i?2 and P2{x) = in the neighborhood of Bi. Multiplying K^^iij} 
by P2 and splitting up in the form ^ = piip + (1 — pi)')p, we have 

l|-f'^0,lV'lloo,So\S2 = \\P2Ko^ii:\\oo,S„\B2 

< C{\\p2KoAPl^A\oo.So\B2 + 11-^^04(1 - Pl)V'l|oc,So\B2) (2-23) 

The first term on the right-hand side of the above inequality contains only an operator 
with a kernel vanishing in a neighborhood of the diagonal x = y, and therefore we 
have 

\\p2Ko,lPlHoo,So\B2<C\mCo{So)- (2.24) 

Since the operator Kq^i mapping C{So) into C"'"(5o) is bounded (see |1|), we find, 
on noting that 1 — pi vanishes in a neighborhood of x* , that 

\\KoA^ - Pimoo,S„\B2 < C\\{1 - Pl)^\Uso\B2 < CMcoiSo)- (2-25) 
From ((2:23)) - ((2:25)1 it follows that 

\\KosH^,So\B2<C\mc,,iSo)- (2.26) 
A similar argument as above gives that 

\\So,i<l>\\oo,So\B2<CU\\co{So)- (2.27) 
Combining (PT^ - d^TTf)) with (P?^ and (P?^ yields 

Moo,s„\B2 < ^^(ll/llco(So) + ll.9llco(So)). (2.28) 

Before proceeding to estimate dv/dv we establish the following estimate in the 
spaces of Holder continuous functions for (V', 0, x? "y^)- 

ll'0llci-°(So\S3) + ll'/'llco-°(So\i33) + llxllc«-(ro) + llv'llc-^iri) 

< C{\\J\\ca(So) + ll.9llco(So) + ll/l|l.a,So\Bi + ||.9l|o,a,So\i3i ): (2.29) 



10 



X. LIU, B. ZHANG 



where is a ball of radius and centered at x* with ri < < We choose 
a function p3 e C^{Sq) such that /03(a;) = for x G So\B2 and P3(x) = 1 in the 
neighborhood of B3. We also choose another function p4 G C^(S'o) such that p4(a:;) = 1 
for X e So\B2 and y04(a;) = in the neighborhood of B3. Splitting U up in the form 



U 











P3(t> 


+ 


(1 


- P3)0 


X 









\ V ) 






1 



■= Up., + C/(l-p3) 



and using Wp^ to denote the matrix W with its first and second rows multiplied by 
P4(a;), it follows from (|2.17p that 



Ro 



^PiU(l-p,). 



(2.30) 

but with two different 



Arguing similarly as in deriving the estimate for ||w||cx;.So\S2 
cutoff functions and replacing pi and p2, we obtain from (|2.30p and (|2.20p that 

||C^P4llo,a = IIV'llco.°(So\-B3) + ll0llco>°(So\-B3) + llxllco.°(ro) + Il<y5||cn.°(ri) 

< C(||i?p,||o,« + \\Ap,Up,\\o^a + ||Ap,[/(i_^3)||o,„) 
<C(||i?p,||o.„ + ||f/||co) 

< C{\\f\\co{So) + ll3llco(So) + ll/lll,c,So\Bi + ll.9llo,a,So\Si), (2.31) 

where || • ||o,q and || • ||co denote the corresponding norms in the product spaces. Now 
it remains to prove the estimate of ll^/iUci, 0(5^^33). Multiplying (|2.13p by P4{x) we 
obtain, on using (|2.30p and noting the fact that the integral operators mapping C"'" 
functions into C^'" functions are bounded, that 

IIV'l|l,a,So\i33 < IIP4V'l|l,a 

< C[\\p4{XoKofl - -ftro,l)V'||l,a + ||P4('5'o,0 " SQs)(t>\\l,a 

+ \\p4{Kofl - M75o,o)x||l,Q + \\p4SaS + 'i'nKo,lS^)ip\\i^a + l|P4/||l,a] 
<C^[||C^IICo + ll(l-P3)t/||o,a+||p4/||l,a] 

< C'lll^^lICo + M\co.''{So\Bs) + U\\c«-''{So\B3) + llxllc«.°(ro) 

+ ll^llco."(ri) + ||/||i,a,So\S3]- (2.32) 



Combining and (lOTjl - fO^ yields the desired estimate (1^:^ . 

. From (|2.12p and the jump relation it is seen 

0,Q,So\S2 



We now estimate 



dv 



that on 5*0, 




dv 






Writing 






dv 




dv 



(2.33) 



Writing ^ ~ p^^ + {1 — p^)ilj andcj) — p^cj) + {1 — p^)4>, we obtain from (|2.33p that 

dv 



< 



0,a,So\B2 



Pi 



dv 



0,a.,Sa 



< C^[ll^llco(So) + U\\co{S„) + 11(1 - P3)V'l|l,a,So 

+ 11(1 - /33)0||o,a,So + ||xllco.»(ro) + ll<^llco."(ri)] 
<C^[IIV^IIco(So) + ll0llco(So) + 11^-111 + \m 

+ llxllco-(ro) + \Mc«■'-{^^)]■ 
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Combining this with (P?^ and (P?^ yields 

< C^(ll/llc„(5o) + ll5llco(5o) + ll/Ill.a,5o\S, + |l.g||o.a,So\Bj.(2.34) 

This completes the proof. □ 

3. The inverse scattering problem. Following the ideas of 12J for transmis- 
sion problems in a homogeneous medium and of [13] for transmission problems in an 
inhomogeneous medium, we prove in this section that the interface 5*0 can be uniquely 
determined by the far field pattern. Combining this with the earlier result in |18j , we 
have in fact proved that both the penetrable interface Sq and the impenetrable obsta- 
cle ^2 with its physical property can be uniquely determined from a knowledge of 
far field pattern. To establish the uniqueness result for the inverse problem, we need 

the following two lemmas, in which = ]R'^\rio so that O2 C and = M'^\J7o for 
some domain VIq with the interface Sq = dflo H dfl and with the domain ^2 C il. 
Lemma 3.1. Suppose the positive numbers fco, fci and Aq are given. For C 

fl, fl2 C let G he the unbounded component of \ (f2 U 17) and let u°°{x,d) = 
u°°{x,d) for all x, d € S with u°°{x,d) being the far field pattern of the scattered 
field u" {x , d) corresponding to the obstacle ^2, the interface Sq and the same incident 
plane wave d). For z ^ G let (u'*, v) be the unique solution of the problem 





Au" + klu" = in f^o \ {z}, 


(3.1) 




Av + klv ^0 in fli, 


(3.2) 




dul dv d<i>o{x,z) 
), -5 = ^ onAo, 


(3.3) 




S§{v) ^OonS*!, 


(3.4) 




lim r{— ikou^) — 0. 

r— *oo OT 


(3.5) 


that is the 


unique solution of the problem — i 


TB with 



0.0, rii, 5*0, Si, SS replaced by fio, Oi, Sq, Si, S^, respectively. Then we have 

u''{x;z)=u%x;z), x eG. (3.6) 




Remark 3.2. By Theorem 12.31 the problem (|3.ip - (|3.5p has a unique solution. 
Proof. By Rellich's lemma [1], the assumption u°°{x, d) = u°°{x, d) for all x, d & 
S implies that 

u^{x,d) =u\x,d), xeG,deS. 

For the far field pattern corresponding to incident point-sources we have by Lemma 
3.3 in |18J that 

$°°(d,z) = $°°(d,z), zeG, 
Thus, Rellich's lemma [4 implies that 

u*(x; z) — u^{x; z), x € G. 
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□ 

Lemma 3.3. Assume that f E L'^{Q,i),h G C{Tq),p G C(ri) and g,ri e C{Sa) 
with r] ^ and rj < on Sq. Then the following problem has a unique solution 

u e c'^{^l)f^C{^nl): 



Au + kfu^ f in r^i, (3.7) 
du 
dv 



— +iriu^ g on Sq, (3.8) 



Furthermore, there exists a constant C > such that 



u ~ h on Fq, (3.9) 
— — \'i\u = p on Fi. (3.10) 



hlloo,n, < Cm\LHn,) + ||5lU,5o + Moo.T, + IbllooTj 



Proof. We first prove the uniqueness result, that la, u ^ Q ii f = g = h ^ p ^ Q. 
With the help of the equation (|3.7p and the boundary conditions p.8p - (|3.10p . we have 



= / {{An + klu)u} dx 
Jni 

/" r iT-^ i2 ,2\ i2i 7 f —du , f _du f _du 

~ / — Vu + Ki m \ dx + j u—ds— / u—ds— / u—ds 



[ {-\\/u\^ + kf\u\'^} dx ~ i [ 7j\u\^diyds + i [ X\ufds. 

Jui Jso JT 



Taking the imaginary part of the above equation, we get u = on some part F of 5*0 
since both rj =/= and 77 < on Sq and A is a nonnegative continuous function. By the 
boundary condition (|3.8p it follows that u = du/dv = on F. Thus, u = in ili by 
Holmgren's uniqueness theorem [H] . 

To solve the problem by means of the integral equation method we introduce the 
volume potential 



{Vf){x) := / <^,(x,y)f{y)dy, 
Jill 



which defines a bounded operator V : L'^{fli) H^{ili) (see Theorem 8.2 in |3]). 
Now look for a solution in the form 

uix) = -{Vf)ix)+ f ^iix,y)[My) + My)]dsiy) 

J So 

^^^^"^^y^ -ij^,{x,y)}{xi{y)+X2{y))ds{y) 



dv{y) 



y) + il Q^^y^ \ (^i(y) + ^2{y))ds{y), x e (3.11) 



with six densities (jji e C{Sq), 4)2 e H2{Sq), xi e C{Tq),^X2 e i?2(Fo), ipi £ C(Fi), 
ip2 € -ff^(Fi) and a real coupling parameter 7 7^ 0. By S we denote the single-layer 
operator (|2.8p in the potential theoretic limit case ki — 0. 
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Then from the jump relations we see that the potential u given by p. lip solve the 
boundary value problem (|3.7p — (|3.10p provided the six densities satisfy the following 
system of integral equations: 

01 - 2{k'q^^ - i'nSo.i)(t>i - Mo(xi + X2) - No{ipi + ip2) = -2g on So, (3.12) 

02 - 2(i^;,i - i7?5o,i)02 = -2 (I; + (Vf) on So, (3.13) 

Xi + 251,1(01 + 02) + MiXi + Ni{ipi + ip2) = 2h on Tq, (3.14) 

X2 + MiX2^2Vf on To, (3.15) 
ipi ~ 2{kI^ + iA5i,i)(0i + 02) - Af2(xi + X2) - N2ipi = -2p on Ti, (3.16) 

^2 - iVi^2 = -2 (^|- + iXj (Vf) on Ti, (3.17) 

where 

Mo = 2{fofi - iiiiKq q + 177^^0,0 + IvSofl), 

No = 2(i^;,i + 77To,i52 + i775o4 - 7?yi^o,i^'): 

A/i = 2(i^l,o-^7^l,o), 
A^i = 2(5l,l+^7i^l,l^2), 

M2 - 2(fi,o - «7^1,o + *Ai?i,o + A7^i,o), 

N2 - 2(^1,1 + ^7fl,lS^2 _ 1^^^2 _^ _ \-iKi^iS^). 

Precisely, we seek a solution {u, 0i, 02, xi, X2, (pi, </'2) G C(r2i) x C(S'o) x i/^ (6*0) x 
C(ro) X H^-{To) X C(ri) X H^{Vi) F to the system of integral equations ([XTT]) - 
(|3.17p . Similarly as in the proof of Theorem 12. 3i it can be shown by using the Riesz- 
Fredholm theorem that this system has unique solution in the space Y and there 
exists a positive constant C > such that 

Il«lloo,ar<^(ll^/lloo,nr + ll5lloo,so + 

+ l|/^l|oo,ro + l|p||co,rJ 
< C(||/|U.(o,) + ||5lloo,So + ll^lloo,ro + lblloo,rj. 

The lemma is thus proved. □ 

We are now in a position to state and prove the main result of this section. 

Theorem 3.4. Suppose the positive numbers ko, ki and Aq fAo 7^ IJ are given. 
Assume that So and So are two penetrable interfaces and £^2 a^ic? £72 are two impene- 
trable obstacles with boundary conditions SS and SS, respectively, for the corresponding 
scattering problem. If the far field patterns of the scattered fields for the same incident 
plane wave u^{x) — e'*^"^ '' coincide at a fixed frequency for all incident direction d ^ S 
and observation direction x d S, then So — So, ^2 = ^^2 ond = . 

Proof. We just need to prove that So = So since the remanning part of the theorem 
then follows from this and Theorem 3.7 in [TB]. Let G be defined as in Lemma [3.11 
Assume that So 7^ 5*0. Then, without loss of generality, we may assume that there 
exists zo e 5'o\f7. Let B2 be a small ball centered at zo such that ^2 nil = 0. Choose 



|; + i,)(V7) 
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h> such that the sequence 



zo + -1^(20) 

J 



.7-1,2,. 



is contained in G D B2, where 1^(2:0) is the outward normal to at zq. Using the 
notations in Lemma [3TT] and letting Vj) and {uj,Vj) be the solutions of (|3.ip - (|3.5p 
with z = Zj. Then, by Lemma [3.1l Uj = Uj := Uj in G. Since zq has a positive distance 

from ^l, we conclude from the well-posedness of the direct scattering problem that 
there exists C > such that 



||Uj||oo,SonB2 + 



9m 1 



< C for all j > 1. 



(3.18) 



oo,SonB2 



Choose a small ball Bi with center zq which is strictly contained in 32- Since 

\\M-,Zj)\\co(So) + ll*o(-,2j)||l,a,So\Si < C", 
li-^(-,2j)llco(So) + \\-Q^i-,Zj)\\o,a.So\Bi < C 

for some positive constant independent of j, we conclude from Lemma |2 .41 that 

dvi 



l|Wjl|oo,So\B2 + 

From this it follows that 



< C for all j > 1. 



',5'o\-B2 



dv 



< G 



for aU j > 1, 
for all j > 1. 



(3.19) 
(3.20) 



00,So\-B2 

The transmission boundary conditions yield 



\\Vj - $o(-,2;j)||oo,SonB2 = 




|oo,SnnB2 










dv 


oo,SonB2 


dv 


00,500-82 



dv 

Combining ([3?20l) and KT2\ yields 

dvj d'^o{-,Zj) 



< G for all 3 > 1, (3.21) 

< G for all j > 1. (3.22) 



An 



< G for all j > 1. 



00, So 



dv dv 

This can be used together with (I3.19P to prove the estimate 

IIAqWj - $o(-,2;i)||co,So < C for all j > 1. 



(3.23) 



(3.24) 



In fact, choose a non-positive function rj e C'^(S'o), ^ and supported in So\B2 
Then Wj := Apfj — $o(', -s^j) solves the following boundary value problem: 

Awj + k\wj — (/cq — fcJ)$o(-,Zj) in fii, 

(^1/7 ' d 

+ iyywj = (— + i'q)[XoVj - $o(-, Zj)] on S'o, 

= -$o(-,2j) on Fo, 



9Wj 



d 

lAwj = -(— +iA)$o(-,^i) on Fi. 
dv 
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Since, by HJIM and (IX^ . / (fcg - fc2)$„(.^ 2^.) ^ L^{n^), h -$o(-,Zj) e C(ro), 
p -(I; + iX)M; e C(ri) and 5 (1; + *^)[Aoi'j - ^oC', z,)] e C(5o), then 
the desired result (|3.24p foUows from Lemma 

Now the triangle inequality together with (|3.2ip and (|3.24p implies that 

ll(Ao - l)^o{-,Zj)\\oo,SonB2 < ||Ao$o(-,Zj) - AoWj lU.SonSa + \\>^oVj - ^o{-, Zj)\\oo,SonB2 

< C. 

This is a contradiction since Aq ^ 1 and \^o{zo, Zj)\ao.SonB2 ^ oo as j ^ cxd. The 
proof is thus complete. □ 

Remark 3.5. Our method can be extended straightforwardly to both the 2D 
case and the case of a multilayered medium, and a similar result can be obtained (that 
is, all the interfaces between the layered media as well as the embedded obstacle can 
be uniquely determined). 
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